Efficient tomography of quantum processes 
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We show with explicit formulas that one can completely identify an unknown quantum process 
with only one weakly entangled state; and identify a quantum optical Gaussian process with either 
one two-mode squeezed state or a few different coherent states. In tomography of a multi-mode 
process, our method reduces the number of different test states exponentially compared with existing 
methods. 
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PACS numbers: 03.65.Wj, 42.50.-p, 42.50.Dv 

Introduction. — One of the basic problems of quantum 
physics is to predict the evolution of a quantum system 
under certain conditions. For an isolated system with 
a known Hamiltonian, the evolution is characterized by 
a unitary operator determined by the Sclrrodinger equa- 
tion. However, the system may interact with its environ- 
ment, and the total Hamiltonian of the system plus the 
environment is in general not completely known. The 
evolution can then be regarded clS cL "black-box pro- 
cess" [l| which maps the input state into an output 
state. An important problem here is how to characterize 
an unknown process by testing the black-box with some 
specific input states, i.e., quantum process tomography 
(QPT) 0. 

Any physical process can be described by a completely 
positive map e. Such a process is fully characterized 
if the evolution of any input state pin is predictable: 
Pout — e(pin)- In general, QPT is very difficult to im- 
plement in high dimensional spaces, and, more challeng- 
ingly, in an infinite dimensional space, such as a Fock 
space. Recently, Ref. Q showed QPT in Fock space, for 
Continuous Variable (CV) states. Two conclusions Q 
are: (i) If the output states of all coherent input states 
are known, then one can predict the output state of any 
input state; (ii) By taking the "photon-number-cut-off 
approximation" , one can then characterize an unknown 
process with a finite number of different input coherent 
states (CSs). 

Here we study QPT using a different approach. Based 
on the idea of isomorphism [3j] , and using the standard Q- 
representation in quantum optics, we show, with explicit 
formulas, that one can complete QPT with either only 
one weakly entangled state for any quantuirr process, or 
only a few CSs for quantum optical Gaussian processes. 
The method described here has several advantages. First, 
it presents explicit formulas without any approximation, 
such as the photon-number-cut-off approximation. Sec- 
ond, it requires only one or a few different states to char- 
acterize a process, rather than all CSs. Third, for multi- 
mode Gaussian process tomography, the number of input 
CSs increases polynomially with the number of modes. 



rather than exponentially. Fourth, it uses Q-functions 
only, which is always well-defined for any state without 
any higher order singularities in the calculation. 

Isomorphism and process tomography with one weakly 
entangled state. — Define 1$^) = X]fc=o 1^^) 
level maximally entangled state in the composite space 
of modes a and b. (Here l^"*") is not normalized, be- 
cause this simplifies the calculations below). Assume 
now that the process e acts on mode b. Using isomor- 
phism 0, if the state = /® e(|<I>^)($^|) is known, we 
shall know £(pin) for any single- mode input state pin on 
mode b. Consider a single-mode input state on mode b, 
\ip{{ck})) = X^fc'^fel^)- Obviously it can be written as 



= tra(ir)(n®/-i*+)(*+i) 



(1) 



and \'ip*)a = J2k Cfc|^)a, is a single-mode state on mode a 
(sometimes we omit the subscript a or b for simplicity). 
The output for any initial state pin = (|V')('*/'l)f, is 

Pout ^ e {tT,[{\r){r\)a® h ■ {\^+){<f+\)ab]} (2) 

Since the partial trace and the map e are taken in differ- 
ent subspaces, their orders can be exchanged. Thus 



Pout=a{r\I^£{\'^^){'i'^\W)a=a{r\Pe\r)a 

= tra(|V*)(^*|®/-Pe). 



(3) 



Equation ([3]) predicts the output state of any input state 
of an unknown process, given p^. However, generating 
the maximum entangled state 1$"*") is technologically dif- 
ficult, especially when s is large. Moreover, for the case 
of CV states in Fock space, s is infinite and the max- 
imum entanglement does not physically exist. There- 
fore, we cannot really test a process with |$+} in Fock 
space. However, we can first test a process with some 
other easy-to-manipulate states, and then deduce the 
state pe. For example, one can test the one-sided map 
/ (g) £ with an arbitrary non-maximally entangled state 
\<Pi{'rk})) = Z]fc=o "rAkk), if ro • ri • • • rs-i ^ 0. Denoting 
the output state as f^{rfc}i we have 



^{r.}=I®e{mrk})mrk})\ 



(4) 
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On the other hand, we know that \4'{rk)) = T{{rk}) ® 
/ !<!>+); and T{{rk}) is a projection operator defined as 
^({^fc}) = Sfc=o Since T®I and the one-sided 
map I ® e commute, then Eq. ([4]) can be written as 

^{r,} = T({rfc}) ®l{l®e [|$+)(<i>+] } f{{ru]) ® I (5) 

which gives rise to 

p,^f-\{rk})®in{.r,}f-\{rk})®I. (6) 

Here f-^{{rk]) is defined as T({rj:^}). If we test the 
one-sided map I ® e with the hmited entangled state 
0({rfc})) and we find that the outcome state is i^{rfc}j 
then, using Eq. ([6]) we can determine the output state 
when the input state is l*!*^). We can then use Eq. ^ to 
predict the output state of any single-mode input state 
on mode h. Explicitly, if the input state is \i^{{ck\)) = 
^^Cfc|fc}, the output state becomes 

p.^ = tra [|^(4/rfe))(V'(4/rfc)| ® I ■ ^{r,}] (7) 

where \il^{cllrk)) ^J2ki4/rk)\k). 

For a state in Fock space, s is infinite and |A;) is a Fock 
state which can be generated by the creation operator 

on the vacuum state |0). We can implement a similar 
technique to the one presented above to characterize an 
unknown quantum optical process with only one weakly- 
entangled state, i.e., a two- mode squeezed state (TMSS). 

Process tomography with one TMSS. — A TMSS is de- 
fined by \x{q)) = Cq exp{qa' b'^)\00) , where Cq = i/l — g^, 
and q is real. The (un-normalized) maximally-entangled 
state here is 1$+) = limq_n exp(ga'''6''')|00) = J2T=o \^^)^ 
where is the creation operator for mode x. 

We define the projection operator T(g) = 
Cq exp[(ln(7)a^a] which has the property: 
f{q) (a, at) f-^{q) = (a/q,gat). The TMSS can 
be written as 

\x{q))^f{q)®I\^+). (8) 

Assume now that the black box process acts only on mode 
b of the bipartite state \x{q))- After the process, we ob- 
tain a two- mode state fig. We now wish to predict the 
evolution of any state under the same process, using the 
information on how the input state Ixil)) changes under 
this map. According to Eq. we have 

nq = f{q) ®I-pe- f{q) ® I (9) 

where pe = I ® e (|<I>+)($+|). Naturally, 

p,^f-\q)®I-nq-f-\q)®I. (10) 

We now also formulate the output state of any single- 
mode input state |'/'({cfc})) = J2k^k\k) of mode b. Ac- 
cording to Eq. ([3]), we obtain the output state 

p^ = tra [\r{{ck/q''})){r{{ck/q''})\ ® I ■ ^q] 

= a{r{{ck/q''}m,\r{{ck/q''}))a. (n) 



More explicit expressions can be obtained by using the 
Q-function. If the single-mode input state on mode b is 
a coherent state |a), the output state then becomes 

Pa = {a*\p,\a*) ^ {a*\f-\q)®I ■nq-f~\q)®I\a*}. 

Note that the state \a*) here is a single-mode coherent 
state on mode a. Using the property of T{q) and the 
definition of CSs, a\a*) — a*\a*), we easily find 

f-\q)®I \a*) ^N'q{a)\a*/q) (12) 

where the factor Nq{a) — exp [— |Q!p(l - l/g^)/2] /cq, 
and \a* / q) is a coherent state on mode a defined by 
a\a* /q) = [a* / q)\a* / q) . Thus, the output state of mode 
b is 

p„ = |A/',(«)|' a{a*lq \^q \a*lq)a ■ (13) 

Assume the Q-function for flq is QQ^{Z*,Z^,ZaTZi,). 
According to its definition, Qci^{Z*,Z^,Za,Zb) — 
{Za, Zi,\D,q\Za, Zh), whcrc \Za,Zf,) is a two-mode coher- 
ent state defined by (a, b)\Za,Zb) = {Za, Zb)\Za, Zb). 
Hence the corresponding density operator is 17^ =: 
Q{a\b\a,b) : , where the normal order notation : . . . : 
indicates that any term inside it is reordered by plac- 
ing the creation operator in the left. For example, 
: aba^'b^'a := a'^b'^a^b. Therefore, using Eq. (|13p and the 
normally-ordered form of flq, we have the following sim- 
ple form for the Q-function 

Qp^iZb*,Zb) = \Nq{a)fQn^{a/q, Zb\ a* /q, Zb) (14) 

of the output state pa- Eqs. ([T31IT4)) are the explicit ex- 
pressions of the output state for the input of any coherent 
state \a). According to Ref. if we know the output 
states for all input CSs, we know the output states of 
all states in Fock space. In our approach, given any in- 
put state I '(/'), we can write it in its linear superposition 
form in the coherent state basis, and then obtain the Q- 
function of its output state by using Eq. (|14|) . These and 
Eqs. (ffl llip . can be summarized as follows: 
Theorem 1: Any process in Fock space is fully char- 
acterized by the bipartite state fig, which is the output 
of the initial TMSS \x{q))i if 9 0- Any process on s- 
dimensional states is characterized by the bipartite state 
^{rk}^ 'v^hich is the output state from the initial bipartite 
state |(?!>({?'fe}), if Tfe 7^ for all fcs. 

Characterizing a Gaussian process by testing the map 
with a few CSs. — One can also choose to test aprocess 
with only single- mode states. As shown in Ref. y , if we 
only use CSs in the test, the tomography of an unknown 
process in Fock space requires tests with all CSs. Though 
this problem can be solved by taking the photon-number- 
cut-off approximation, in a quantum-optical process as- 
sociated with intense light, one still needs a huge number 
of different CSs for the test. Here wc show that the most 
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important process in quantum optics, the Gaussian pro- 
cess, can be exactly characterized with only a few CSs in 
the test. 

A Gaussian process maps Gaussian states into Gaus- 
sian states. Therefore the Q-function of the operator 
must be Gaussian: 

Qp^{z:,z;,Za,Zi,)=cMco+L+L^+S+S^+So), (15) 

Zn \ ^ I ^ n \ I Z„ 



where L = Q 

So = {z:,z;)Y 



Za 



Zb 

-^aa -^ab \ ^j^^ Y I ^ab 

Xba Xbb J ' \ Yba Ybb 

Before testing the map, ah these are unknowns. The 
normally-ordered form of the density operator is : 
Qp^{a\h\a,h) : . The output state from any single- 
mode input coherent state \u) (on mode b) is 



S = ^{Za,Zb)X 

g = (r^Tb), X ^ X 



Zb 

T 



Pu = tra [(|w*)(u*|)a ®I-Pe 

Its Q-function is 



(16) 



QpuiZb^Zb) = Qp^{u,Z^,u*,Zb) 
= exp{cu + Lu + Ll + R + R^ +Ro), (17) 

where Lu = {Tb + u*Xab + uYab)Zb, R = ZbXbbZb/2, 
Rq ~ Z^YbbZb, and Cu is determined by cq, Fq, Xaa and 
Yaa- Exphcitly, 

C„ = Co + Re {2TaU* + U*XaaU* + uYaaU*) (18) 

The quadratic functional terms, {R, R\ Rq) on the ex- 
ponent in Eq. (jl7p are independent of u; these terms must 
be the same for the output states from any initial CSs. 
Therefore, these can be known by testing the map with 
one coherent state. Thus, we do not need to consider 
these terms below. Now suppose that we test the pro- 
cess with six different CSs, \ai), and i = 1, • • • 6. Assume 
also that the detected Q-function of the output states is 

Qp^^ {Z;,Zb) = exp(c, + D, + dI+R + R^ + Ro) (19) 

where Dt — diZb is the detected (hence known) lin- 
ear term. According to Eq. {17} . the Q-function of the 
output state from the initial state \ai) of mode b must 
be Qp^.{Z^,Zb) = Qp^{ai,Zl,ai*,Zb). Therefore, we 
can derive self-consistent equations by using the detected 
data from p^^ and setting u ~ Ui in Eq. (|17p : 



Li Ca^ Ci 



(20) 



where Li, Cq. are just Lu, c«, respectively, after setting 
u ~ ai in Eqs. p7m8|) : Di and q are known from tests. 
Explicitly, Li — (Tb + Oii* Xab + OiiYab)Zb- The first part 
of Eq. (pO)) causes: 



1 ai* ai 

where A' = | 1 a2* C(2 \ , d = {di, d2, d^)^ . There 
1 as* "3 

are three unknowns (Tb, Xab, and Yab) with three equa- 
tions now. We find 



(Tb, Xab, YabV = K~^d . 



(22) 



If the Gaussian process is known to be trace-preserving, 
then Eq. ([^ completes the tomography: up to a nu- 
merical factor, we can deduce all the output states of the 
other input CSs, \ai), for i = 4, 5, 6. The term Ci can be 
fixed through normalization, which is determined by the 
quadratic and linear functional terms on the exponent of 
the Q-functions. Knowing these {q}, one can construct 
Pe completely, as shown below. For any map, c; can be 
known from tests with la,). We then have 



'-^ ' {^0, ^a, Fq, Xaa, X aa, ^aa) 



1 „,* 



J 



1 0^2 Ct2 

1 al as 

1 a| a4 

1 al as 

V 1 "6 "6 



2 "3 
2 4 
2 "5 
2 "6 



for the second part of Eq. (f20| . Thus 



C2 
C3 
C4 



(23) 



(24) 



(co. Fa, Fq, Xaa, Xaa, ^a) J C 



Theorem 2: Given K and J defined by Eqs. 
then the QPT of any single-mode Gaussian process in 
Fock space can be performed with six input CSs, when 
deti^ ^ and det J ^ 0. The QPT of any trace- 
preserving single-mode Gaussian process in Fock space 
can be executed with three input CSs, when det K ^ 0. 

For example, one can simply choose ai = 0, a2 = 1, 
aa = i, = —1, = —i, and ag = 1 + i. One finds 



Co 


= ci, Fb = di 






Xab 


= [-{l + i)di 


+ d2 + ids] /2 




Yab 


= [-{l-i)di 


+ d2- ids] /2 




Fa 


= (C2 + iC3 - 


C4 - ic^) /4 




Xaa 


= [2ici + (l- 


2i)c2 - (1 + 22)C3 + C4 - C5 - 


h 2ic6] /4 


Y 


= -Ci + (C2 4 


- C3 + C4 + C5)/4 


(25) 



where {di} and {ci} are defined in Eq p9)) . 

An example. — As a check of our conclusion, we calcu- 
late the output state of a beam-splitter (BS) process as 
shown in Fig. 1. The BS has input modes b and c and 
output modes b' and c'. Regarding this as a black-box 
process, the only input is mode 6 and the only output is 
mode b'. We set mode c to be vacuum. The BS trans- 
forms the creation operators of modes b, c by: 



K ■ {Tb, Xab, Yab) 



(21) 



Ubs (&^ ct) Ug^s = {b\ ct) M 



BS 



(26) 
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Beam-splitter 
b 




Input of the 
Gaussian Map 



Output of the 

Gaussian Map 



FIG. 1: Gaussian Map constructed by a beam-splitter 



where Mbs 



If we test such a pro- 



cos 6 sin 
— sin 9 cos 9 

cess with a coherent state la^), we shall find Pq. = 
I cos 9ai) (cos 9ai \ . Comparing this with Eq. (|19p , we have 
di = a* cos 9 and ci = —\aiC0s9\'^. Using Eqs. 
we find 



bb 



-1, Xab - 



cose, Ya, 



bb 



COS 

= Co 







(27) 



Therefore =: exp(a^5t cos 6* — atacos^6' — 6^5 -t- 
ab cos 9) : . With this we can predict the output state 
of any input state, for example the displaced squeezed 
state |C(r, Z)) = exp (-§6+2 + r^2) ^^^t _ |o), 
where r is real. According to our Eq. ([3]), 
Pi = tra[(|?(r,Z*))(e(r,Z*)|)a^4-Pe]. As a resuh, 
Qp,{Z;, Zb) = Cexp(Hi -1^2+^3 + 7^4), where C is 
the normalization factor, and Tii = |.Zf,p(tanh^ r sin^ 9 — 
% = (^6^ + zf)tanhrcos2 0/(2(7), - 
ZbCose(Z* - Ztanhrsin2 6')/(gcoshr), Hi = 
Zj;'cos0(Z - Z*tanhrsin^0)/(gcoshr), and 
g = 1 — tanh^ r sin^ 9. This is same with the re- 
sult from direct calculations using Eq. (pS)) . 

Multi-mode extension. — Multi-mode Gaussian QPT 
has many important applications. For example, it applies 
to a complex linear optical circuit with BSs, squeezers, 
homodyne detections, linear losses, Gaussian noises and 
so on. Consider now a Gaussian process acting on a k- 
mode input state (on mode 61, 62, • • • ^fe), with outcome 
also a fc-mode state. Even though other methods 0] can 
also be extended to the multi-mode case, the number of 
input states required there increases exponentially with 
the number of modes fc, because the number of ket-bra 
operators in Fock space increases exponen- 

tially with k. As shown below, the number of input states 
in our method increases polynomially. 

To apply isomorphism Q , we consider k pairs of maxi- 
mally entangled states, each on modes ai,&i; 02, &2,-'" 
ak,bk. Explicitly, |$+) = |0+)2 • • • |</)+)fc. Here 



= lim„ 



exp ^qalbf^ |00) indicates a maximally- 



entangled state on modes at, bi. Subspaces a and b each 
are now fc-mode. Any state lip) in subspace &, can still 
be written in the form of Eq. ([l}, with the new defi- 
nitions for |-0) and \^'^)- Using Eq. (jlOp . it is obvi- 
ous that the output state of these /c-pairs-TMSS fully 
characterize the process. A fc-mode QPT can also be 
tested with fc-mode CSs, if the process in Gaussian. The 
main Eqs. (|22l [24]) still hold after redefining the nota- 



tions there. First, Fa, Fb, u, ai, di, Za,Zi, are now 
/c-mode vectors. For example, \ai) = {an, ai2, ■ ■ ■ ctik), 
Zb = [Zbi, Zb2, - ■ ■ Zbk), di = {dii,di2, - ■ ■ diu), and so 
on. Following Eq. ()15p . Xxy is now a kxk matrix, for 
X = X or Y with x = a,b y = a,b. We still apply 
Eqs. (EllMl) to calculate {Ffc, Xab, Yab} and {F^, X^a, 
Yaa}i respectively, but keep in mind that the matrices 
K, J and symbols d, c are now redefined. There are 
{2k + l)k unknowns in (Fs, Xab, Yab)- We need 2fc + 1 
different CSs of fc-mode to fix these unknowns. Matrix K 
is now {2k + 1) x {2k + 1), since each ai here is a fc-mode 
row vector. Here d is a (2fc + 1) x fc matrix as d'^ = 

(df, dTi---c^2A;+i)' "^itli '^i = {dii,di2, - ■ -dik). Simi- 
larly, J is now iVx and N — (fc-M)(2fc-|-l), since af and 
|aip here are row vectors of af = {En, Ei2, ■ ■ ■ Eik) and 
\ai\^ = {En, Ei2, - ■ ■ Eik), and each element of Ei„i (or 

Ei„i) is a vector with (fc — m -I- 1) modes (or fc modes), as 
,2 



Eini 
Eini — 



{a 
{a 



1 ^ini^i ^m-\-l ; 



i7n ' 

imO^il: ^tm^i2i 



a. 



ak,a\) and 
). Obvi- 



ously, c is a column vector with N elements. Therefore 
we conclude with this: 

Corollary 1: Any fc-mode map e in Fock space is char- 
acterized by the output state of fc-pair-TMSS under one- 
sided map I (E) s. Any fc-mode Gaussian QPT can be 
performed with (fc -I- l)(2fc -|- 1) different CSs of fc-mode; 
or with 2fc -I- 1 different CSs of fc-mode if the process is 
trace- preserving. 

In summary, we have presented explicit formulas quan- 
tum process characterization with only one weakly en- 
tangled state, as well as the tomography of a quantum 
optical Gaussian process with a few different coherent 
states. These results have been extended to multi-mode 
quantum optical process and the number of test states 
required increases only polynomially with the number of 
modes. 
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